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Abstract — The economic dispatch problem is considered for 
unbalanced three-phase power distribution networks entailing 
both non-deferrable and elastic loads, and distributed generation 
(DG) units. The objective is to minimize the costs of power 
drawn from the main grid and supplied by the DG units over a 
given time horizon, while meeting the overall load demand and 
effecting voltage regulation. Similar to optimal power flow coun- 
terparts for balanced systems, the resultant optimization problem 
is nonconvex. Nevertheless, a semidefinite programming (SDP) 
relaxation technique is advocated to obtain a (relaxed) convex 
problem solvable in polynomial-time complexity. To promote a 
reliable yet efficient feeder operation, SDP-compliant constraints 
on line and neutral current magnitudes are accommodated in the 
optimization formulated, along with constraints on the power 
factor at the substation and at nodes equipped with capacitor 
banks. Tests on the IEEE 13-node radial feeder demonstrate the 
ability of the proposed method to attain the globally optimal 
solution of the original nonconvex problem. 

Index Terms — Unbalanced distribution systems, economic dis- 
patch, power factor, voltage regulation, elastic loads. 



I. Introduction 

The advent of distributed energy resources, along with 
the rapid proliferation of controllable loads such as, e.g., 
plug-in hybrid electric vehicles (PHEVs), call for innovative 
energy management methodologies to ensure highly efficient 
operation of distribution networks, effect voltage regulation, 
and facilitate emergency response [1]. Toward these goals, 
variants of the optimal power flow (OPF) problem have been 
devised with the objective of optimizing the power supplied 
by distributed generation (DG) units as well as by the utility 
at the substation, subject to electrical network constraints on 
powers and voltages, and the expected load profile [2]. 

These approaches however, are deemed challenging because 
they require solving nonconvex problems. Non-convexity 
stems from the nonlinear relationship between voltages and 
the apparent powers demanded at the loads. Furthermore, 
the high resistance-reactance ratio in conventional distribution 
lines severely challenges the convergence of Newton-Raphson 
iterations, which have been traditionally employed for solving 
nonconvex OPF problems in transmission networks [3], [4]. 
This has motivated the adoption of forward/backward sweep- 
ing methods [5], which enable computationally-efficient load 
flow analysis, but are not suitable for optimization purposes, 
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fuzzy dynamic programming [6], particle swarm optimiza- 
tion [7], sequential quadratic optimization [8], and steepest 
descent-based methods [9]. However, these approaches gen- 
erally return sub-optimal load flow solutions, and may be 
computationally cumbersome. To alleviate these concerns, the 
semidefinite programming (SDP) reformulation of [10] and 
[11] was recently extended to balanced distribution networks 
in [12], and conditions ensuring global optimality of the 
obtained solution were derived. 

Three-phase distribution feeders however, are inherently 
unbalanced because i) a large number of unequal single- 
phase loads must be served, and ii) non-equilateral conduc- 
tor spacings of three-phase line segments are involved [13]. 
As a consequence, optimization approaches can not rely on 
single-phase equivalent models (as in, e.g. [12], [9]). For the 
unbalanced setup, an OPF framework was proposed in [14], 
where commercial solvers of nonlinear programs were used, 
and in [15], where Newton methods where utilized in con- 
junction with OpenDSS load flow solvers. A model based on 
sequence components was adopted in [16], and the Newton- 
Raphson algorithm was used. However, since these methods 
are inherently related to gradient descent solvers of nonconvex 
programs, they inherit the limitations of being sensitive to 
initialization, and do not guarantee global optimality of their 
solutions. 

The main contribution of the present paper consists in 
permeating the benefits of SDP relaxation techniques to the 
economic dispatch problem for unbalanced three-phase power 
distribution systems. This powerful optimization tool not only 
offers the potential of finding the globally optimal solution of 
the original nonconvex problem in polynomial-time complex- 
ity [17], but also facilitates the introduction of thermal and 
quality-of-power constraints without exacerbating the problem 
complexity. The focus here is on the case where the costs of 
power provided by the utility company and supplied by the DG 
units are known in advance over a given time horizon. Then, 
the goal is to minimize the overall energy cost so that both 
non-deferrable and elastic load demands are met, and the node 
voltages stay within prescribed limits. Furthermore, constraints 
on line and neutral current magnitudes, as well as on the power 
factor at substation and nodes equipped with capacitor banks 
are accommodated in the optimization problem in order to 
improve reliability and efficiency of the distribution feeder. 

Notation: Upper (lower) boldface letters will be used 
for matrices (column vectors); for transposition; (•)* 

complex-conjugate; and, complex-conjugate transposi- 

tion; 5R{ } denotes the real part, and 5{ } the imaginary part; 
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j = represents the imaginary unit. Tr(-) denotes the 

matrix trace; rank(-) the matrix rank; o the Hadamard product; 
and, | • | denotes the magnitude of a number or the cardinality 
of a set. Given a vector v and matrix V, [v]-p denotes a IT^I x 1 
sub-vector containing the entries of v indexed by the set V, 
and [Vj-p-^-p,, the |Pi| x \T J 2 \ sub-matrix with row and column 
indexes described by V\ and Vi . Finally, Omxn and 1m xn 
denotes MxN matrices with all zeroes and ones, respectively. 

II. Modeling and problem formulation 

Consider a radial distribution feeder comprising N nodes 
collected in the set TV := {1,...,7V}, and overhead or 
underground lines represented by the set of edges £ := 
{(m,n)} C (TV U {0}) x (TV U {0}), where the additional 
node corresponds to the point of common coupling (PCC). 
The feeder operation is to be optimized over a given time 
interval X := {1, 2, . . . , T}, where each time slot can represent 
e.g., ten minutes or fifteen minutes, one hour, etc, depending 
on the specific short-, medium-, or long-range scheduling 
horizon [14]. 

The backbone of the feeder generally consists of three-phase 
lines, with two- and single-phase connections at times present 
on laterals and sub-laterals. Let V m n C {a, b, c} and V n C 
{a, b, c} denote the phase of line (to, n) G £ and node n G TV, 
respectively. Further, let V„ t G C and if t G C be the complex 
line-to-ground voltage at node n G TV and time slot t of phase 
<fi G V n , and the current injected at the same node, phase, and 
time. As usual, the voltages vo,t := [Vg* t , V^ b t , Vj^] 7 " at the 
PCC are assumed to be known [13]. 

Lines (to, n) G £ are modeled as 7r-equivalent components, 
and the |P mn | x |P mra | phase impedance and shunt admittance 
matrices are denoted as Z mn G c' 7 '" 1 "^' 7 ''""' and G 
(Q|P m „|x|;p m „| ; respectively. If four-wire grounded wye lines 
or lines with multi-grounded neutrals are present, matrices 
Z m „ and Ymn can be obtained from the higher-dimensional 
"primitive" matrices via Kron reduction [13, Ch. 6]. Using 
the 7r-equivalent model, it follows from Kirchhoff's current 
law that the current lf % can be expressed as [13] 



- V 

mgAV, 



r(s) 



[V„, t ]-p„ 



V m ,t TV 



(1) 

where TV„ C TV is the set of nodes linked to n through 
a transmission line, and v n t G C' 73 "' denotes the column 
vector collecting the voltages at node n and time slot t. Three- 
or single-phase transformers (if any) are modeled as series 
components with transmission parameters that depend on the 
connection type [13, Ch. 8], [14]. 

Per phase <fi G V n an d node n G TV, the following two 
classes of loads are considered. 

• A base non-deferrable load with active and reactive 
powers demanded at time t denoted by p£ t and t , 
respectively. 

• A set T)^ of controllable (elastic) loads, each with a 
prescribed energy requirement E% n to be completed over 
a given interval zj n := [4, n JdJ Q z < with 4,n 
representing the starting time, and f$ n the termination 



slot; that is, J2 t£xtn Pt n ,At = K* with P tn,t the 
amount of active power supplied to the controllable load 
d G V% at time slot t, and A t > the duration of the 
time slot. 

An example of controllable load is PHEVs, whose charg- 
ing process can be shifted from hours with high price of 
electricity [18], and high load conditions of the distribution 
network [8], [19], to off-peak hours. In this case, users specify 
the time when PHEVs will be plugged in, and the time by 
which the charging has to be completed [8], [18]. 

In distribution feeders, capacitor banks are mounted at some 
selected nodes to provide reactive power support, aid voltage 
regulation, and correct the load power factor (PF). As usual, 
capacitors can be modeled as wye or delta loads with constant 
susceptance [14], [13, Ch. 9]. Therefore, with yf, denoting 
the susceptance of a capacitor connected at node n and phase 
<j>, the reactive power Qq n t provided by the capacitor at time t 
is given by Qq t = yf, n \V^ t \ 2 . To satisfy the load demand, 
DG units such as, e.g. diesel generators and fuel cells can 
be employed to complement the power drawn from the main 
distribution grid. Then, suppose that S DG units are located 
at nodes S c TV, and let t and Q^ G t denote the active 
and reactive powers supplied by unit s G S. 

The focus is on the case where the costs of power provided 
by the utility company P$ t := 3t{V£ t (I$j)*}, (j> = a,b,c, 
and supplied or consumed by the DG units are determined 
in advance for the period X. Let {«o,t} denote the former, 
and {c s t } the latter. Then, the goal is to minimize the 
overall cost of power purchased from the main grid and 
generated within the feeder (economic dispatch), so that the 
total load demand is met, and the node voltages stay within 
prescribed limits; that is, the following problem is to be solved, 
where V« := {{ijj^.t, {< t , < t , ^,n,*> <„>#,»,t} 
and Vd '■= {{Pdn t}v0,n,t} collect the optimization variables: 

(PI) min£ L,t £ Po t + E C M E P L,t) (2a) 

V V d tex \ 4>ev sgS <j>ev s J 

s x.v* t (itr = p*-p LiS > 



+ jQc,s,t - 3Qtn,V V t G X, <t> G V., S G S (2b) 

V <t> (T<t> \* _ _p0 _ p<l> _ AC)^ 

v n,t\ 1 n,t) — r L,n,t r d,n,t J^L,n,t 

d£T>t 

+ jyc,JV* t \ 2 , VteX l( f>er n ,neAf\S (2c) 
E P d,n,At - E+ n , V d G X>t 4> 6 V n , n G N (2d) 

< P|„ )t < P d m ™ VteI,deVtneM (2e) 
v™r < \V* tt \ < vtex^e P n , n G TV (2f) 
PZ?.<P&,., t <P%£, VteX^eV^seS (2g) 

Qg% < Qc,s,t <Qg% V t G X, 4> G V n , n G TV (2h) 

where P^ p g^^Qg%Qg% capture physical and opera- 
tional constraints of the DG units, and VJf lln and V" laK are 
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given minimum and maximum utilization and service voltages. 
Finally, P™jf represents a possible cap for Pf n . If capacitor 
banks are not present, d2cb should be modified to V„ t (I„ { )* = 
- p Ln t ~ JQl n v Recall that the voltages at the PCC are 
assumed known. However, if needed, this assumption can be 
relaxed, and (PI) can be appropriately re-stated. 

Unfortunately, (PI) is a nonlinear nonconvex problem due 
to the load flow equations (f2bb-(f2cb as well as the voltage con- 
straints ( f2fb - In the next section, an equivalent reformulation 
of (PI) will be derived, and its solution will be tackled by 
employing an SDP relaxation technique. 

III. Relaxed semi-definite programming 

Consider first a distribution feeder with only three-phase 
lines and nodes; that is, \V n \ = 3 for all n E Af, and \Vnt\ = 
3 for all lines (n,l) € £. Let Y G cHN+i)x3(N+i) be a 
symmetric matrix defined as [cf. (Q]i] 



[Y] 



Er 



V n .V„ 



_7-l 

inn ? 
^3x3: 



±Y 



inn 



if m = ii 

if (to, n) G E 
otherwise 



and define the Z(N+\) x 1 vectors v t := [vq t , . . . , t ]' r and 
i t := [i£ t , . . . , il >t ] T , with i n , t := [Il t ,I^t, Il t ] T . Then, © 
can be re- written in vector-matrix form as i t = Yv t . 

Since the PCC voltages {vo.t} are known, re-write the 
vector of complex voltages as v< = ao t o x t , with x t := 



L 3 ) v l,i> • ■ • i V N,t\ and a 0,« — L v O,t> J -|P l |) • • • , *-\V N \} 

for all tel. Then, consider expressing the active and reactive 
powers injected at each node at time t, as well as the voltage 
magnitudes as linear functions of the outer-product matrix 
X t := x t x^. To this end, define the following admittance- 
related matrix per node n and phase (f> 



IT 



Yt :=ei(ei) r Y 

where := [0^, . . . , 0^,, e^ r , Oj^, . . 
and {e^}0 g { a c } denotes the canonical basis of 
for future use the Hermitian matrices 



(3) 
Denote 



' U |-P«| 



■ P.n.t 



iD«(Y r t + (Y ? t)«)D , t 



Q,n,t 



V.n.t 



2 

■i-T) n (Y 4 
TyH p*(p0)T n 



(4a) 

(4b) 
(4c) 



with Do t := diag(ao,t). Using a linear model in X t 
(and therefore in V t := v t v^) is established in the following 
lemma (see also [20] and [11]). 

Lemma 1: Apparent powers and voltage magnitudes are 
linearly related with {X t } as [cf. ©J 



Tr(*£, n ,tX t ) = P*„,t - Pt,n, 



r Q,n,t^t) — Qa.n,t 

Tr(*t, n ,tXt) = |< t | 2 



Q 



V p 

/ , 1 d,s,t 

-^, n Tr(* 



(5a) 



v.n^t) (5b) 
(5c) 



with P£„ >t = Qt, n ,t = for n G and y* „ = if 

capacitor banks are not present at node n. 



Proof. See the Appendix. □ 
Using (f5at — (T5cT>. problem (PI) is equivalently reformulated 
as follows: 



(P2) min V 
{x f} ,v d ^ 



E 



Ko,t 7 . Tr(* P0 1 



Xt) 



(6a) 



S .tTr(#* X t ) + P£ 



/ . x d,s,t 



0. 



Tr(*S,«.* X * 



VteX,ct>eV n , VneN\S (6b) 

2t,„,t-^,„Tr(*t,„Xt) = 0, 

Vtel,<j)£ T n , Vn G A/"\S (6c) 
P^<Tr(# PiSit Xt)+P^ Sit + 

V tel,<t>eV s ,VseS 



T}<P s- pmax 



(6d) 



Minn 



rmin\2 



VteI,(t>eV s ,VseS (6e) 

(C)' < Tr(#t,„, t X f ) < (VT ax ) 2 ,V^ Vn e Af (6f) 



E 



„ t A t = \/d&Vt<p&V n ,n£N (6g) 



0<P|„ it <P d % Vtel,deVtneAf (6h) 
rank(Xt) = 1, Vtel (6i) 
Xt >r 0. [X t ] Po ^ = l 3x3) V t G 2. (6j) 

Unfortunately, (P2) is still nonconvex because of the rank-1 
constraint on the positive semi-definite matrices {X t }. Nev- 
ertheless, (P2) is amenable to the SDP relaxation technique, 
which amounts to dropping the rank constraints, thus relaxing 
nonconvex problems to SDP ones; see e.g., the tutorial [17], 
and the works in [20] and [11], where this technique is 
employed for power system state estimation and OPF for 
power transmission systems, respectively. Leveraging the SDP 
relaxation technique here too, it is possible to obtain the 
following convex relaxation of (P2): 



(P3) 



mm 



{x ( }.v d 



E%ETr($t, , t X t ) 



tez 4>ev a 

+ EE c ^E Tr (*p,,tX t ) 

tez ses 4>ev s 

s.t. Xt t o, [XtK, Po 



(7a) 



l 3x3 , and (|6b| - (|6hJ» . 

Clearly, if all the optimal matrices {X° pt } of (P3) have 
rank 1, then the variables {X° pt },V^ pt represent a globally 
optimal solution also for (P2). Further, since (PI) and (P2) 
are equivalent, there exist 2\I\ vectors {x° pt } and {v° pt }, with 
X° pt = x° pt x° ptW and v° pt := ao,t o x° pt , for all t G I, such 
that the optimal objective functions of (PI) and (P2) coincide. 
This is formally summarized next. 

Proposition 1: Let {X° pt }, V^ pt be the optimal solution of 
(P3), and assume that rank(X° pt ) = 1, for all t G X, Then, a 
globally optimal solution of (PI) is given by V^ pt , the vectors 
of complex line-to-ground voltages 



opt 



Vtel 



(8) 
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where Ai,t G M + is the unique non-zero eigenvalue of X° pt 
and Ui.t the corresponding eigenvector, and the supplied active 
and reactive powers 



popt 



Pq% = Tr(*S, M v? Pt vr W ) + Qt s , t ,V S £SU {0}. (10) 



" pis (9) 



The upshot of the proposed formulation is that the globally 
optimal solution of (P2) (and hence (PI)) can be obtained via 
standard interior-point solvers, in polynomial-time complexity; 
see, for example, the complexity bounds for SDP reported 
in [21, Ch. 4] and [17]. This is in contrast with gradi- 
ent descent-based solvers for nonconvex programs, sequen- 
tial quadratic programming, and particle swarm optimization, 
which in general do not guarantee global optimality of the 
obtained solutions, face challenges pertaining to sensitivity of 
the initial guess, convergence, and complexity which increases 
with the number of iterations. Notice also, that matrices 
{$P„ ( , *q „ ti *v>i m ' e vel T sparse. This property can be 
leveraged to substantially reduce the computational burden 
of interior-point solvers; for instance, the so-called "chordal" 
structure of matrices {X f } can be effectively exploited, as 
advocated in [22]. 

Since (P3) is a relaxed version of (P2), matrices X° pt 
could have rank greater than 1. In this case, rank reduction 
techniques can be employed to find a feasible rank-1 ap- 
proximation of X° pt (see [17] and references therein). The 
resultant solution is feasible for (P2), but generally subopti- 
mal [17]. Notably, when balanced distribution networks are 
considered, [12] established conditions on the voltage angles 
and the reactive power injections under which rank-1 matrices 
are always obtained provided the non-relaxed problem is 
feasible. Derivation of similar conditions in the present context 
constitutes an interesting future research direction, that will 
naturally complement the result in [12]. 

A. Feeders with two- and single-phase lines 

For feeders with two- and single-phase laterals and sub- 
laterals, the dimensions of matrix Y have to be adjusted to 
12n=o \Pn \ x J2n=o \Pn\> an( ^ ^ s entr i es have to be as follows: 

i) matrix — occupies the \V mn \ x \V mn \ off-diagonal 
block corresponding to line (m, n) £ £; and, 

ii) the \V n \ x \V n \ diagonal block corresponding to node 
n e TV U {0} is obtained as 



[Y] 



E 



- 1 - ran 



(ii) 



vW if 

1 ran 11 



where Z mn = Z mn and 
V n = Vmn, otherwise [Z mn ] Vnm „v nm = z ™« and 
[^mn]v n \P nm ,V n \P nm = (Ymn is computed likewise). 
Re-defining the J2n=o \P™\ x 1 vectors collecting voltages 
and currents as v' t := [v£ t , [vi,*]^ , . . . , [v Nyt ] VN ] T and 



Lii,*]p i: 



, [ijv.tl-Pjy] 7 "' respectively, ([TJ can be 
re-written again in vector-matrix form as \' t = Yv' t , for all 
t € I, and the procedure (IUi-© can be readily followed. 



IV. Feasible voltage profile 

To effect voltage regulation, and avoid abrupt voltage drops, 
constraints on the minimum and maximum utilization and 
service voltages were imposed in (PI). Constraints (f2fb how- 
ever, may challenge the feasibility of (PI), since it may not 
be possible to meet the minimum (maximum) utilization and 
service voltage requirements when feeders are heavily stressed 
and DG units supply a substantial amount of power. It is thus 
of prime importance to perform preemptive analysis of the 
feasible voltage profile in order to unveil possible infeasibility 
of (PI) and, in case, facilitate corrective actions. To this end, 
consider solving the following optimization problem 

(PA) min (I-W)EE E (Kt\*-\V? i3 ' 



E 




c s ,t 



E^ 



G,s,t 



(12a) 



i>£V s 



with wy € (0, 1) denoting a weighting coefficient, and |V^ f | 
the prescribed voltage magnitude of the feeder (e.g., | V„ f | = 1 
p.u.). Although constraints < f2fb are not present in (P4), the 
first term in ( I12ab promotes regulation by penalizing voltage 
magnitudes that deviate from the nominal ones. 

Similar to (PI), problem (fT2l is nonconvex. However, by 
exploiting again Lemma [T] along with the SDP relaxation 
technique, the following convex problem is obtained 



(P5) 



min (1 — wv) > 
{x t }.v d ' ^ 



t,n,(fi 



+ Wv J2*o,t E Tr (*p,o, t x t ) 

tex 4>ev a 
+ w v E E c ^ E Tr (*L,t X *) d3a) 



tex ses 



s.t. 



~ UL n,t 

Tr(*t.„X t ) - \V, 



ref|2 



Tr(*£, n X t ) 
-1 



■ef|2 



-< 



(13b) 



X t h 0, [X t ]p , Po = 1 



3X3: 



and (|6b) - ® 



where constraint (1 1 3bb is enforced for all nodes n, per phase 
(f>, and time slot t. Notice that by using ( f5cT > the first term 
in (I12ab becomes quadratic in {X t }. To bypass this hurdle, 
the non-negative real variables {a^ t } are introduced to upper 
bound each term (Tr(*^ n X t ) - jv^ 2 ) 2 , and the Schur's 
complement is subsequently employed to obtain ( 113bl i. If for 
a given wy, all matrices {X t } have rank 1, then the optimal 
solution of (P5) is also a globally optimal solution of (P4) (see 
Proposition [TJ. 

Clearly, if the voltages {V^'° vt } obtained from (P5) satisfy 

(Kt) 2 < \Vh° pt \ 2 < (V™r) 2 , for all t,n,<t>, then it is 
possible to proceed with the solution of the economic dispatch 
problem (P3). On the other hand, if some of the voltage 
magnitudes largely deviate from |V^ f |, corrective actions have 
to be taken; these include, for example, switching the taps 
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of controllable capacitor banks, or curtailing portion(s) of the 
loads. 

V. Thermal and Quality- of- power constraints 
A. Thermal constraints 

High current magnitudes on the lines can have detrimental 
effects on both efficiency and reliability of the distribution 
network. From an economical perspective, an additional (real) 
power has to be drawn from the main grid, or, supplied 
by the DG units in order to compensate for the increased 
power dissipated on the distribution lines. On the other hand, 
conductors may overheat if stressed by high currents over a 
prolonged time interval, and may eventually fail. This, in turn, 
would trigger an outage event, with consequent interruption of 
the power delivery in portions of the feeder. 

To alleviate these concerns, it is of interest to constrain 
either the power dispelled on the conductors, or, the magnitude 
of currents flowing through the distribution lines [12], which 
amounts to adding one of the following constraints in (PI): 

l^mn.tl — Iran (14) 
P mn '■= |- r mn| 2 ^{[ Z mn]{0},{0}} < P~mn (I 5 ) 

where I^ nn t and P^ n denote the current flowing on the phase 
cj> of line (m, n) 6 £ , and the active power lost on the same 
line and phase, respectively. 

An SDP-consistent re-formulation of (fl4li-([T"5ll has to be de- 
rived in order to accommodate the aforementioned constraints 
in (P2) and (P3). To this end, let i mn .t := [{I mn t\V denote 
the \P mn \ x 1 vector collecting the complex currents flowing 
through line (m, n) G £, which is related to the line-to-ground 
voltages v„,t and v n ,t as (cf. (l}) 

imn.t = Z~J, ([V m , t ]-p mn - [V n ,t]v m j ■ (16) 

Notice that since is generally not diagonal [23], ( [TBI 

captures also current components arising from mutual induc- 
tive reactances and capacitive coefficients. Define the |P mn | x 
J2n=o com pl ex matrix 

B mn := [0|-p mn | xJ2 m =o l^nl' Z m«' ' • ' 

ip m „ixi:^;L + i i - Pni' z ™«°i^m»ixE^=„ + i \Vtd ( - 17 - ) 

where Z™ n is a \P mn \ x \P m \ matrix with elements 
fc]^„?„ = Z- x „ and [Z™ n ] Vmn ,v m \v mn = 0; likewise, 
ZJ^„ has dimensions |P mn | x \V n \, and its entries are filled as 

KJp m „,P™ = -^mn and fej?™,?,^ = °' Thus ' 
building on ( TToT l. an SDP-compliant re-formulation of (Tflli- 
(TT3T > is provided in the following lemma. 

Lemma 2: Consider the Hermitian matrix 

*/,mn,i := D o!t B mr l e J tr l ( e J tr l ) rB " l nDo,t (18) 

where {& mn } <t>£V m „ denotes the canonical basis of JRl Pmn L 
Then, constraints ([T4l-([T5Tl can be expressed linearly in the 
outer-product X t as 

Tr{<f>t mn t X t } < I%£ (19) 
5R{Tr{et„(etj r Z mn }}Tr{*f m „ t XJ < P™. (20) 



Proof. See the Appendix. □ 

In distribution feeders, line outage events maybe triggered 
by overheating effects on the neutral cable(s), especially those 
experiencing highly unbalanced load conditions. Towards de- 
riving constraints on the magnitude of neutral current(s), 
let Vmn denote the set of grounded neutral cables that are 
present on the line (m, n) <G £, and T„ m the IT 5 ™! x |"P mn | 
neutral transformation matrix obtained from the primitive 
impedance matrix of the distribution line via Kron reduc- 
tion [13, Sec. 4.1]. For example, the neutral transformation 
matrix of a four-wire grounded wye segment has dimensions 
1x3, while its dimensions increase to 3 x 3 for an underground 
wye line with three neutral conductors. Thus, the neutral 
currents i^,t : = [l£l,u ■ ■ •> I mn,h T are linearly related to 
the line currents i m n,t as [13, Sec. 4.1] 

imn.t ~ ^-mni-mn.t • (21) 

It readily follows from (I2TI 1 and the result of Lemma |2] that 
the magnitude of the current on the neutral cables can be 
constrained in the SDP problem (P3) as 

Tr{S$m (t X t } < 4r2< max , V <p G P$£ (22) 

with 

*W T) n Ti n T H ( P (v) ^T T R n c23 , i 

*/,mn,t '~ L, 0,t D mn 1 mn t! mnl ,! iimi mn mn *-> ,t 

where, as usual, {e™]} is the canonical basis of M.^™^. 

There is an increasing concern on the thermal effects arising 
from harmonic currents in the neutral cable(s). In this case, 
constraints similar to (1221 can be imposed on a per-harmonic 
basis (see, e.g. [9]). 

B. Constraints on the power factor 

The PF has been increasingly recognized as one of the 
principal measures of efficiency and reliability of power dis- 
tribution networks [1], [24], [9]. High PF translates to lower 
generation and transmission costs, and enhanced protection of 
transmission lines from overheating (hence, higher resilience 
to line outages). Constraining the PF at the PCC is tantamount 
to limiting the reactive power exchanged with the main power 
grid. This, in turn, has two well-appreciated merits: i) it 
alleviates the power losses experienced along the backbone 
of the feeder [24]; and, ii) it limits the current drawn at 
the PCC, and therefore facilitates coexistence of multiple 
feeders on the same distribution line or substation without 
requiring components such as, e.g. conductors, transformers, 
and switchgear of increased size. 

Unfortunately, the definition of PF for an unbalanced 
polyphase system is not unique [25]. In this paper, a per- 
phase definition is adopted in order to limit the reactive power 
exchanged at the PCC on each phase. Intuitively, polyphase 
variants [25] may induce high discrepancies between the 
amount of reactive power exchanged per phase, but with a 
"good" polyphase PF nevertheless. Let t/q t G [0, 1] denote the 
minimum PF required at the PCC on the phase <fi G Vq at time 



IEEE TRANSACTIONS ON POWER SYSTEMS (SUBMITTED) 



6 



slot t e I. Now consider adding the following constraint to 
(PI): 

-l 



<t (l 



0.1 



\I+ 
\ 1 0,t 



> 



(24) 



where voltages {|Vq t |} are known [13]. Notice that DG units 
complement the power supplied by the utility, and are usually 
not sufficient to satisfy the load demand on their own. Under 
this premise, an SDP-consistent reformulation of (l24l can be 
readily obtained, as summarized next. 

Lemma 3: Provided the power supplied by the DG units 
does not exceed the total load demand at the feeder, d24l i is 
equivalently expressed as a linear function of X as 



< t Tr(*t,„X t )-Tr(** jn X t )>0 
< t Tr(*^X t )+Tr($^„X f )>0. 



(25) 

□ 



Additional charges are generally applied to residential and 
industrial loads with a poor PR In the presence of highly 
inductive loads, capacitor banks are usually employed to 
balance reactive demand, and thus maintain the PF as close as 
possible to 1 [13]. Recall that yf, denotes the susceptance 
of a capacitor connected at node n and phase <j>, and the 
provided reactive power amounts to Qf, n = Vc n \^n,t\ 2 - 
With Q^ L n t > denoting the reactive power demanded by 
an inductive load, a minimum per-phase PF 77^ t £ [0,1] is 
imposed as [cf. d24l i1 



P 



L,n,t 



(Pi 



L,n,t) 



(Q 



L,n,t 



!C,n,t 



>< t V^G?„ (26) 



where Pf 



L t is given. Clearly, |V^ t | can be re-expressed as a 
linear function of X t using ( TScT i. and ( f26b can be reformulated 
to obtain the following SDP-compliant form. 

Lemma 4: Using d5c| l, constraint (|26l l is equivalent to 
the following linear matrix inequality (with Qq n t = 



P 



P 



L,n,t 



Q 



L.n.t 



L.n.t 



-1 



■ L,n,t 



'C,n,t 



-< 



(27) 
□ 



Proof. See the Appendix. 

Controllable capacitor banks can be accounted for by asso- 
ciating an integer variable with each of the capacitor switches. 
To tackle the resultant mixed integer nonlinear problem, ex- 
haustive search over the switches can be performed [14]. In 
this case, (P3) is solved for each switch configuration. 

VI. Numerical Tests 

The proposed optimization framework for unbalanced three- 
phase systems is tested on the IEEE 13-node test feeder shown 
in Fig.Q] Compared to the original scheme in [23], DG units 
are placed at nodes 1 and 10. Specifically, single-phase DG 
units supply a maximum real power of 300 kW and 500 kW, 
respectively, and they operate at a unitary PF. Capacitor banks 

















KM 


si 


t 


3 ; 



{a, V'} 
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fa,»,c} 



(«■'"■} ^ 



Fig. 1. Modified IEEE 13-bus test feeder. 

with rated reactive power of 200 kVAr and 100 kVAr are 
present at nodes 5 and 8, respectively. Line impedance and 
shunt admittance matrices are computed based on the dataset 
in [23]. To solve (P3) (and (P5) for a preemptive voltage profile 
analysis), the MATLAB-based optimization modeling package 
CVX [26] is used, along with the interior-point based solver 
SeDuMi [27]. 

The time horizon is 24h, and slots of lh are considered; 
that is, X = {1AM, 2 AM,..., 11PM, 12 AM}. The loads 
specified in [23] are assumed to be the peak demands of 
the day, and the "spring mid-week" load profiles reported 
in [28] are used to generate {_Pf in t , Q\ „ t }t&x- Specifically, 
the "commercial load profile" in [28, Sec. 1.1] is used for node 
9, whereas the "residential load profile" is applied to all the 
remaining nodes; a Gaussian random variable with mean 1 and 
standard deviation 0.1 is used to insinuate a perturbation on 
the profiles on a per-node and a per-time basis. The resulting 
aggregate real loads per phase are depicted in Fig. [2] (similar 
trends are obtained for the reactive loads, but are not reported 
due to space limitations). 

Ten controllable loads are present at node 5', and are 
allocated as follows: 2 on phase a, 4 on phase b, and 4 on 
phase c. The energy requirement is 11 kWh, and the cap -P™ x t 
is set to 4 kW, so as to resemble the demands of 10 PHEVs [8]. 
Customers are assumed to plug-in the PHEVs at 6 PM, and the 
charging has to be completed by 6 AM. Two additional elastic 
loads are present at node 9, and have to be satisfied between 8 
AM and 4 PM. In this case, the energy requirement per each 
load is 30 kWh, and no cap is present for {Pf g t }- 

To model the price of the power purchased from the main 
distribution grid {ko,*}, one-day ahead locational marginal 
prices (LMPs) available in the Midwest Independent Trans- 
mission System Operator (MISO) [29] database are utilized. 
Specifically, the LMPs for the Minneapolis area on June 7th, 
2012 are utilized throughout this section, and are reported in 
Fig. |2] On the other hand, the cost incurred by the use of the 
DG units is kept constant over time, and it is set to 30 $/MW. 

A minimum PF of T)q t = 0.8 is required at the PCC, and the 
limits V r f m = 0.95 p.u. and V,? lax = 1.05 p.u. are imposed 
to enforce voltage regulation. Finally, a balanced flat voltage 
profile is assumed at the PCC, with |Vjf t | = 1.02 p.u., and 



ZV b t = 120 c 



Constraints 

on the line currents are not considered, since this datum is not 



and AV§ t = -120 c 



IEEE TRANSACTIONS ON POWER SYSTEMS (SUBMITTED) 



7 





Aggregate 


v.;:.; 


3hase a 










■ - * - 

2 r "* " 


Aggregate demand, phase c 
_UP [t/MW] 










o - 








a* . ' * 


a--;' 




»-A ..... 








h -o- - 


J* ~ 


3 □ C 


^J- '*«"-■.> „~ * 

-□- -a- « -j s 4 


it 

6- " 


-0- -0- <T r 
J- -n. _^ 










''' 


•..•mi 


,* 










2 - 












0' 

2 AM 2 


kM 4 


M 6 AM 8 AM 10 


m 12 


PM 2 


M 4 PM 6 PM 8 PM 10 


PM 12 

















Aggregate 


demand, phase a 




















Aggregate 
Aggregate 


demand, phase b 
demand, phase c 


- 
































































































AM 2 


M 4 AM 6 


M 8 AM 


10 


12 


PM 2 


M 4 


M 6 


M 8 


M 10 


PM 12 



Fig. 2. Aggregate real load profile [MW], and prices {^o,t} [$/MW]. 
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Fig. 4. Aggregate per-phase elastic demand [kWh]. 
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Fig. 3. Real power supplied [MW]. 

available in [23]. 

Before proceeding, it is worth mentioning that the rank of 
the matrices {X° pt } was always 1 in the experiments reported 
in this section. Therefore, the globally optimal solutions of 
(PI) were always attained. This illustrates clearly the merits 
of the proposed formulation. 

Fig. |3] depicts the active power supplied by the DG units, 
and drawn from the main distribution grid. As expected, the 
DG units are heavily utilized from 8 AM to 11 PM, which is 
the interval where the price of power purchased from the main 
distribution grid is higher that 30 $/MW. This, in turn, has the 
benefit of reducing the overall demand of the feeder during 
the peak hours (peak shaving). Notice however, that the DG 
units are not utilized at the maximum extent because of the 
constraint on the PF at the PCC [30], as it will be shown later 
on. The optimal allocation of the elastic energy demands is 
shown in Fig. [4] It is shown that the PHEVs connected at node 
5' are charged from 2 AM to 5 AM, interval where both the 
LMPs and the non-deferrable loads are the lowest. A similar 
behavior is noticed for the elastic demands at node 9; in fact, 
they are entirely satisfied in the time slot [8 AM, 9 AM], which 
is the slot with the lowest LMP in the interval [8 AM, 4 PM]. 

Finally, Fig.[5]portrays the trajectories of the PF at the PCC. 
It can be seen that the lower bound on the PF is tightly met 
when the DG units are active. In fact, as they supply real 
power to a lagging power system, a reduction of the PF is 
inevitably experienced at the PCC. This can be further noticed 
from the dotted (orange) trajectories, which correspond to the 
case where (P3) is solved without the constraints on the PF. In 
this case, the majority of the real power is supplied by the DG 
units, thus entailing a significant drop of the PF at the PCC. 
The case where the DG units can operate at a variable PF, 



Fig. 5. PF at the substation. 

from 0.5 to 1, is also considered. In this case, the DG units 
can supply a sufficiently amount of reactive power, so that the 
PF at the substation can be kept close to the unity most of the 
times. This suggests that DG units can be effectively utilized 
for providing reactive support [30], although an appropriate 
modeling of the cost incurred in this case is required. 

VII. Concluding remarks 

The paper considered the economic dispatch problem for 
unbalanced three-phase power distribution networks, where 
the costs of power drawn from the main grid and supplied 
by the DG units over a given time horizon was minimized, 
while meeting the overall load demand and effecting voltage 
regulation. Is spite of the inherent non-convexity of the for- 
mulated problem, the SDP relaxation technique was advocated 
to obtain a (relaxed) convex problem. As corroborated by 
numerical tests, the main merit of the proposed approach 
consists in offering the potential of finding the globally optimal 
solution of the original nonconvex economic dispatch problem. 

Appendix 

Proof of Lemma Q] To prove (Bat , notice first that the in- 
jected apparent power at node n, phase <f> and time t is given by 
VU< t r = = (vM(et) T it) n - Next, noticing 

that v f = an t o x t — Do tx t and using i f = Yv t , it follows 
that (-vMH) T U)" = (x«D«ef(el)^YD , f x f )« = 
(x*D* Y^D , t x t )« = x«D« (Y*)«D 0lt x t , which can 
be equivalently rewritten as Tt(Dq t Y^Dn,tXt). Thus, the 
injected real and reactive powers can be obtained by us- 
ing, respectively, the real and imaginary parts of (YJJ) . 
Finally, Oct can be readily established by noticing that 
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Kt? = v?efM) T v t = x«X>«e*(ejJ)7-D , t x t = 
Tr(D«e£(e^D , t X t ). 

Proof of Lemma [2] From ( fTSI l. and using the defini- 
tions of Z™„ and Z" m , it follows that i m „ jt i^„ t = 

^mn v m,t y m ,t £J mn ^ ^mn* n,t v n.t^mn ZJ mn v m,t» n>i ii mn 

Z"„v n , t < t Z™« = B mn v t v*B*„. Thus, |J*J 2 is 
given by |J&J 2 = (etJ r B m „D , t x t x«D« B«„e*„ = 
Tr{B m „D , t X t D« B«„et n (etJ r } = Tr{<t>? m „ t X t }. 

Proof of Lemma [4] After standard manipulations, d5c| i can 
be re-written as (P* n , t + - y^Tr^^X,)) 2 < 

(P| n t r]~ j) 2 , which is quadratic in X t . Then, ( |27] i is readily 
obtained by using Schur's complement. 
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